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The topic of this paper is to study the two way time-transfer problem between a ground based 
station and a low orbit Earth's satellite, in the aim of an application to the T2L2 experiment. 
The sudy is driven in a fully relativistic framework. Because of the rapid increase in clock's preci- 
sion/measurements, the first term beyond the Earth's potential monopolar term is explicitly taken 
into account. Explicit formulae, for both the distance and offset problems (definitions in the text) 
are proposed for the relevant applications. 
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I. INTRODUCTION 

The idea of optical space-based time transfer had first 
been proposed to ESA in the seventies through the 
project which was named LASSO p] (Laser Synchroniza- 
tion from stationary Orbit). Even though it was quite 
difficult to get some data with LASSO, the mission was 
a great success with some time transfers obtained within 
a time stability in the range of 100 ps. In 1994 OCA 
proposed to build a new generation of optical transfer 
based on some new technology: it was the beginning of 
the T2L2 project. As compared to LASSO, the objec- 
tive was to improve the performances by at least 2 or- 
ders of magnitude. After several proposals on the Mir 
Space Station, ISS, GIOVE, and Myriade, T2L2 was fi- 
nally accepted in 2005 as a passenger instrument on the 
Altimetry Jason-2 satellite ([2], [3]). 

The full exploitation of the T2L2 experiment requires a 
relativistic description of the Earth gravitational field for 
both the definition of the proper timescales relevant for 
the (ground-based and on board) clocks and the photons' 
trajectories involved by laser links. 

The aim of the current paper is to propose explicit 
useful formulae dealing with two-way experiments, like 
T2L2. We derive explicit formulae giving the ratio Vq / Uq 
of the ground-based frequencies at reception and emission 
respectively as a function of the required quantities (po- 
sitions, velocities,...) at a ground-based time. We also 
give explicit expressions for the time dependance of (1) 
the ground proper time needed for a photon to go to the 
satellite and to come back and (2) the offset between the 
on board and the ground based clocks. All the quantities 
are explicitely written as a function of the time attached 
to a unique event (emission or reception). To anticipate 
the performance increase of future clocks, these formulae 
take the J2 term in the Earth gravitational field effect on 
the photons' orbits into account. 

Blanchet et al [3] derive explicit formulae in both the 
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one-way and the two-way problems. However, in the two- 
way case, the ratio Vq/vq of the ground-based frequen- 
cies at reception and emission respectively and other use- 
ful expressions are not explicitly given as a function of 
the required quantities at a unique time. The current pa- 
per proposes such formulae, which can be directly used in 
practical applications. Besides, the photon's orbit is de- 
termined in [J] using the monopolar Earth gravitational 
term only. (A more complete gravitational modelization 
is used in [5], but for the link between proper and coor- 
dinate timescales only.) In the present paper, we include 
the J2 gravitational contribution to light propagation, 
which will be useful as soon as the 10 -14 s level of preci- 
sion is required, as it should be in the future ACES exper- 
iment [5], for instance. Using the Synge world function 
approach [B], Linet and Teyssandier [7] have obtained 
explicit expressions of laser link in the one-way case for 
general axisymetric gravitational field, but no formulae 
in the two-way case is given. This problem has been re- 
considered by Le Poncin-Lafitte et al [8], but they give 
explicit formulae in the spherical case and for the one-way 
problem. Let us remark that, from a theoretical point of 
view, the applicability of the Synge world function ap- 
proach can be limited to regions of spatial extension of 
the order of 2R Sun {R Sun / 'R Sc hw) ~ 2000 A.U. ~ 0.01 pc 
in the worse case (light rays grazzing the Sun and consid- 
ering points at typically ~ 1000 A.U.). The last remark 
is obviously not relevant for solar system applications, 
but is relevant for interstellar scales. 



II. EXPERIMENTAL DESCRIPTION 

A. Principe 

The experiment principle is issued from satellite laser 
ranging techniques. T2L2 permits to synchronize remote 
ground clocks and compare their frequency stabilities 
with a performance never reached before. T2L2 allows 
synchronization of a ground and space clock and mea- 
surement of the stability of remote ground clocks over 
continental distances, itself having a time stability in the 



2 



Clock (USO) 



Photo detection system & 
tagging unit 




^return with the following time equation: 



Laser station 



FIG. 1: Ground to space time transfer with T2L2 
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FIG. 2: Global synoptic of the T2L2 instrument on Jason-2 
together with the LRA and the DORIS USO 



range of 1 ps over 1000 s. 

The principle is based on light pulses propagation be- 
tween ground laser stations and a satellite equiped with 
a photo detection system and a time-tagging unit (fig. 
fl]). The space instrument uses a Laser Ranging Array 
(LRA) which is also used for some laser ranging pur- 
poses and the Doris ultra-stable oscillator (USO) funda- 
mentally used by the DORIS positioning system (fig. [2J. 
The USO is linked to the time tagging unit. It is the 
onboard reference clock of the project. 

A laser station at ground emits short asynchronous 
laser pulses towards the satellite. The LRA returns a 
fraction of the received photons back to the station. To 
describe roughly the main idea of the principle of the 
measurements, the station records the start (t s tart) and 
return (i rotU rn) times of each light pulse. The T2L2 pay- 
load records the arrival time (iboard) in the temporal ref- 
erence frame of the onboard oscillator. These data are 
regularly downloaded every 2 hours to the ground via a 
regular microwave communication link. 

For a given light pulse emitted from ground, the time 
offset #gs between ground clock G and space clock S 
is deduced from the measurement triplets £ s tart, Aboard, 



(^start "t" ^rcturn)/2 Aboard ~T" 



(1) 



where e is a corrective term, including atmospheric effects 
for instance. The time transfer between several laser sta- 
tions at ground (A, B, C, ...) is then deduced from the 
difference between each ground to space time transfer. 

However, a formulae like ([l]) is based on some inad- 
equate physical concepts from the start. Obviously, it 
involves only a rough notion of time, leaving in the dark 
the fact that different proper times have to be considered 
in the problem, as it is well known from the relativistic 
gravity theory. In fact, ([T]) is correct as long as the new- 
tonian description of (space-)timc correctly describes the 
physics. However, at the level of precision reached by ex- 
periments from years ago, it is well known it is no longer 
the case. Of course, it is possible to put the so-called 
" relativistic corrections" into the e corrective term in ([!]) . 
However, it is better to reformulate the problem directly 
in the relativistic framework, this way ensuring that all 
the "relativistic effects" will be consistently taken into 
account at the considered order. 

Depending on the distance between the laser ground 
stations (roughly 6000 km for Jason2) T2L2 can be oper- 
ated in a common view mode or in a non-common view 
mode. In common view configuration, with two laser 
ranging stations A and B firing towards the satellite si- 
multaneously, the noise of the on-board oscillator has to 
be considered over time interval equal to the delay be- 
tween consecutive laser pulses. In such a way, the noise 
coming from the onboard oscillator can be considered 
negligible in the global error budget. In a non-common 
view mode, the temporal information is carried by the 
satellite local oscillator over the distance separating the 
two ground stations visibility and the USOs noise be- 
comes an important part. 



B. T2L2 on JASON2 

Jason- 2 is a French- American follow-on mission to 
Jason- 1 and Topex/Poseidon. Conducted by NASA and 
CNES, its goal is to study the internal structure and dy- 
namics of ocean currents mainly by radar altimetry. For 
the needs of precise determination of the satellite orbit, 
three independent positioning systems are also included: 
a Doris transponder, a GPS receiver and the LRA (Laser 
Ranging Array) retroreflector. The satellite was placed in 
a 1,336 km orbit with 66° inclination by a Delta launcher. 
The time interval between two passes varies from 2 to 14 
hours with an average duration of about 1000 s. The 
T2L2 instrumentation contains two main subsystems: of 
a photo detection module and an event timer (fig. |3j. 
The photo detection unit is made with 2 avalanche photo 
diodes. One is working in a special non-linear Geiger 
mode for precise chronomctry. The other works in lin- 
ear gain mode in order to trigger the whole detection 
chain and to measure the received optical energy and the 
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FIG. 3: Synoptic of the T2L2 space instrument 



continuous noise flux originating in the earth albedo. To 
minimize the false detection rate, the detection threshold 
may be adjusted either by remote control or automati- 
cally as a function of the continuous flux measurement. 
Some elements of the photo detection unit are located 
outside the main Jason-2 payload on a boom which sup- 
port also the LRA. These elements are oriented toward 
the earth surface. The rest of the detection is, together 
with event timer, inside an unique package located inside 
the Jason-2 payload. All this equipment has a mass of 
10 kg and a power consumption of 42 W. 



C. T2L2 Objectives 

The objectives of the T2L2 experiment on Jason-2 are 
threefold: 

• Validation of optical time transfer, including the 
validation of the experiment, its time stability and 
accuracy. It should further allow to demonstrate 
one-way laser ranging 

• Scientific applications concerning time and fre- 
quency metrology allowing the calibration of ra- 
diofrcquency time transfer (GPS and Two- Way), 
fundamental physics with the measurement of light 
speed anisotropy and alpha fine structure constant, 
Earth observation and very long baseline interfer- 
ometry (VLBI). 

• Characterization of the on-board Doris oscilla- 
tor, especially above the South Atlantic Anomaly 
(SAA). 



D. First results 

T2L2 relies on the international laser ranging network 
that permanently contributes to the International Earth 



Reference System (terrestrial frame and Earth's rota- 
tion), and to the traking of several scientific space mis- 
sions. Since 1998, this network has been organised as 
an International Service (named ILRS), as it was the 
case for the GPS and VLBI geodetic space techniques. 
Among 35 Satellite Laser Ranging (SLR) stations, 17 
of them provide full rate ranging data on Jason-2 that 
are downloaded to the European Data Center (EDC in 
Muechen, Germany) regularly; the data files are up- 
loaded by the T2L2 Scientific Mission Center, located 
in Grasse, France. Among these stations, 10 of them use 
the new data format that permits to extract the start 
time of laser pulses at the ps level, whereas the others 
use the current format. This data format was fitted for 
laser ranging only and use a maximum of 12 digits to 
record the start epoch of the laser pulses (corresponding 
to 100 ns). 

Several steps are necessary to proceed together ground 
SLR and on board T2L2 data sets (the T2L2 data files 
are uploaded from CNES, Toulouse to the Mission Center 
every day with a time delay of 1-2 days): 

• selection of data relatively to the SLR stations and 
satellite passes, 

• extraction of measurement triplets from both data 
sets, 

• estimation of instrumental corrections, 

• determination of the time of light travel between 
the SLR station and the T2L2 space instrument, 

• computation of the ground to space time transfer. 



1. Data selection 

The first step of the data treatment consists in deter- 
mining the absolute frequency offset and the phase delay 
between the space clock and the UTC time scale, along 
the time. The goal is to compute the on board date of 
each detected optical event (i.e., arrival time of emitted 
laser pulses) in a time scale that is identical to UTC at 
the picoseconde level. Because the second step of the 
data treatment (Triplets extraction) , is based on a direct 
comparison of the dates of T2L2 data and ground SLR 
data. 

The T2L2 space instrument being connected to the 
GPS of the Jason-2 plateform, we use Pulse Per Second 
(PPS) that are permanently emitted by this system at 
il^isec and then recorded by T2L2 in the space clock 
time scale. That permits to establish a phase link and 
an absolute frequency offset between the space clock and 
the UTC-GPS time scale, respectively to 0.3 /isec and to 
2-3. 10" 11 Hz (over 5000 sec). 
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2. Triplets extraction 

The second step of the data treatment consists in rec- 
ognizing from all the events recorded by T2L2 those cor- 
responding to emitted laser pulses. Because the events 
coming from the SLR stations are blended together, it 
is necessary to find from which station a given event is 
coming up. On the other hand, due to the combination of 
several parameters such as (i) the low energy level used by 
current SLR systems (around 50-100 mJ and even lower), 
(ii) the energy dispersion due to atmospheric effects, (hi) 
and the aperture size of the T2L2 optical module rela- 
tively to the diameter of the laser beam, all the emitted 
laser pulses don't have a corresponding on board event. 
As a statistical result, T2L2 is detecting more events from 
SLR stations that use a higher energy level for their laser 
pulses (there is a factor 2 to 3 between the stations). 

3. Instrumental corrections 



35 ps depending of each station. By taking into account 
the cinematics between the satellite and the ground sta- 
tion during a satellite pass, it is possible to compute a 
synthetic time of flight by averaging the row data locally 
over an integration duration of a few tens of seconds. 
The noise is thus decreased by a factor 1/y/N where N 
is the number of row data. If that computation is done 
properly from a data set to another it permits also to 
reduce the slight but significant differences between the 
laser stations. 

The noise of the time of light travel being reduced to 
4-5 ps inside the time transfer, the time stability can 
be represented by a computed time variance between 
roughly 1 s and 600 s. As an example, the time sta- 
bility of the Wettzell's Hydrogen Maser compared to the 
T2L2's DORIS quartz oscillator is 40 ps at 1 s and 7 ps 
at 30 s. For time integration greater than 30 s this mea- 
surement is limited by the DORIS time stability which 
is 5 ps at 30 s and 10 ps at 100 s. 



Several instrumental corrections must be estimated 
both for the ground and space segments. At ground, the 
accuracy of the range is obtained by a calibration pro- 
cess which uses a calibration target located at a known 
position by a classical survey at the mm level. At the 
satellite, several parameters have to be taken into ac- 
count as (i) the geometrical delay between the reference 
points of the LRA and the T2L2 optical module, (ii) the 
time walk compensation of the photo detection module 
which is sensitive to the number of received photons, (hi) 
the event timer calibration based on an internal measure- 
ment process. 

All these corrections are computed at the ps level, 
using the available mechanical, optical and electronical 
measurements that were determined for all Jason-2 sci- 
entific equipments before its launch. In addition, the 
stellar sensors of the satellite are providing the attitude 
informations needed to compute the geometrical effects. 

4- Ground to space time transfer first results 

The ground to space time transfers that have been se- 
lected use a set of SLR data from very performant ground 
stations that, in addition, are equipped by a Hydrogen 
Maser or a Cesium clock as a permanent system of time. 
Thus, the present analysis takes into account data sets 
from Wettzell (Germany), Matera (Italy), Changchun 
(China), and Grasse (France) stations. 

Nevertheless, if all these selected SLR stations provide 
data with an optimal precision (range and start times), 
a noise is introduced in the time transfer by the time in- 
tervals measured at ground. It mainly comes from the 
quality of the laser beam (i.e., the pulse width), the 
ground detection module (a photo detector), the opti- 
cal behaviour of the on board LRA, and the atmosphere 
effects (two times); the current budget error is of 20 to 



5. Improvements 

Now, we can routinely compute time transfers between 
different ground clocks and the space DORIS oscillator. 
Nevertheless, some work is still required to improve the 
instrumental model (physical behaviour of the hardware) 
to the level of 1-2 ps at 30 s. As examples of this model, 
(1) the T2L2 photo detection module is sensitive to the 
number of incoming photons that should be properly cor- 
rected, and (2) there is still a data filtering to develop in 
order to avoid Solar events that are detected by T2L2 
just 1-3 ns before the real pulse event. 

But the current estimation of the T2L2 performances 
seem to be in a good accordance with the previous speci- 
fications of the overall project thus justifying the present 
relativistics equations to be applied. 



III. THE SPACE-TIME MODEL USED FOR 
TIME TRANSFER 

In GR, photons follow null geodesies of the space-time 
metric g^. We write the metric under the form 

ds 2 = gooc 2 dt 2 + 2goicdtdx l + gijdx 1 dx 3 , (2) 

where ct = x°. In the IAU2000 recommendations, terms 
in c -4 are neglected in while they are taken into ac- 
count in both goo and goi (with no contribution in goi) 
[3]. Since we are interested in light motion, dx ~ cdt, con- 
sidering c -4 terms in goo (and goi) would be meaningless. 
Hence, the metric to be considered writes 

ds 2 = (-1 + ^jf\ ^d^-^cdtdx'+^l + <%cteW 

(3) 
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where [9] 
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a and a 1 being the gravitational mass and mass current 
respectively. ||u|| is a short notation for V u k u k . Hence, 
w corresponds to the potential in Newtonian gravity. Let 
us write w under the form 
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where we is the Earth potential, written as the sum of 
the monopolar term WE, m , the oblateness term We,j 2 , 
and the other terms we, UP : when developed in spherical 
harmonics, and w ext the external potentials, caused by 
the presence of the Moon and the Sun essentially. The 
orders of magnitude of these different terms are given by 
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where Me and Re are Earth's mass and radius, and U 
the Newtonian potential for a remote body of mass M, 
at distance L. It results that, for satellites at low alti- 
tude (say < 1000 km), the term WE, m induces corrections 
which can reach ~ 10 -11 s for the time transfer. The ef- 
fect of the J2 term can then be of the order of 10~ 14 s. 
The following terms in the Earth potential development 
results into corrections of order 10~ 16 s. For both the 
Sun and the Moon, w ex t is of order 10~ 19 . 

The time-space term in the metric (dragging effect) is 
of order 
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where p is the Earth's density, and V the field of matter 
velocity inside the Earth. Assuming a spherical rigidly 
rotating Earth, with uniform density and angular veloc- 
ity (jj, the expression ^ leads to, for a satellite located 
in the equatorial plane at altitude h 
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(7) 



for a low altitude satellite. The resulting correction to the 
photon flying time is then expected to be of the order of 
Hr 17 s. 

In a foreseeable future, time transfer precision of order 
10~ 15 s is expected. Hence, for low altitude satelites, one 
has to consider the J 2 term in Earth potential while all 
the other terms may be discarded. 



IV. APPLICATION 



A. The photon orbit 



From the previous section, we will consider in this sec- 
tion the time transfer problem in a space-time described 
by the metric 



>h 2 = f -1 + J) c 2 dt 2 + 7 ^ <),,, /,•<,/,-/ (8) 



where 7 is a PPN parameter including possible extensions 
of GR, like scalar-tensor (ST) theories [TU]- The potential 
w is given by 



w + J 2 w 



with 
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and 



w 



GM E R 2 
2r 3 



(9) 
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R being a parameter of the order of the Earth radius. 
Since we are interested in photon propagation, one has 
to solve the geodesic equation 



1 



dX 
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with 



k a 



dx c 



(11) 

(A being an affine parameter). The wave quadri- vector k 
satisfies the null condition 



k a k p = 



(12) 



expressing the motion occurs on light cones. The way 
chosen to solve equations |TT] ) and (12), up to first order 

in w, is the way chosen in |llj . in which only the lu-term 
was considered in the potential. The aim is to derive 
the trajectory of the photon x % (x°), obtained after the 
elimination of the affine parameter. (This elimination 
can be made both after the integration of the equations 
or at the level of the differential equation ( jTT| ) itself.) 

Up to 0(c~ 3 ) terms, the result reads 
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where, using the short notation m = GMe/c 2 
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suitable for a satellite and one can replace W by w in 
this case in practical calculations). One has 
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where is the position at £ = (the integration con- 
stants Z % being adjusted this way [H]). The three num- 
bers rf are also integration constants satisfying the nor- 
malization condition 



VLU 1 = 1. 



In the first order terms, one has 
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and an analogous expression for dr B . Since t B = ^a + £, 
where £ is the time transfer (to be determined later), this 
leads to 



dr B 
dr A 
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The indexes recall that the corresponding quantities 
have to be computed at the first photon emission time. 
The function £ (£ e ) as a function of the emission time £ e 
(= £a) is implicitely obtained from the equation defining 
the interception of the photon by the station B 

x ph (£, n k , x\ (£ e )) = x B (t, x% (£ e ) , v k B (£ e ) , ...) (22) 

where one has explicitely written the dependence of 
the station B in initial (i.e. at the emission time) or- 
is ex- 



bital parameters. The function x ph (t 



(*e)) 
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plicitely given by (131, (|14j) up to first order. The func- 



(£ e ) (£eJ7.-) can be obtained from the 
geodesic equation if the station B is a satellite while it is 
given by an Earth rotation model in the case of a ground 
based station. Let us write these functions as 



f = 4 - n^zg. (19) 
Let us remark that = K 2 and n l £ l = 0. 



B. Time transfer and useful formulae 

1. The time transfer 

Let us now consider two photons emitted from the sta- 
tion A and received by the station B. Let drA and dr B be 
the proper time intervals between, respectively, the two 
emissions and the two receptions, and dtA, dt B the corre- 

be the 



AO' ^BO 



sponding coordinate time intervals. Let x^ 
positions and v A0 , v B0 the coordinate velocities of A and 
B at the coordinate time corresponding to the emission 
of the first photon by A. 

Since both A and B have velocities corresponding to or- 
bital free motions (for a satellite) or smaller (for a ground 
station), the velocities of both are at best of the order of 
y/GMs/r, in such a way that the g oi term in (Jij) con- 
tributes as a c~ 4 term in each station proper time. The 
ground station's potential has to be modelized more pre- 
cisely than in ([£]) |3] . Hence we generally denote it by W 
and not by w which is its approximation ([9| (w is only 



x ph (t,n k ,x k A (t e )) = x A (t e ) + n l ct 



+ x ph (*> nk t x A (*e)) 

+ ( % 2)l ph (t,n k ,x k A (t e )) (23) 



and 
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(24) 



x B (£ e ) + tv B (£ e ) + \t 2 a B (£ e ) + \t 3 b B (£ e ) 
2 



where 



(3/2)' 



ph 



corresponds to the (unknown) 3/2 order 
part of the photon orbit (contribution of the <?oi term). 
The coefficients a B (£ e ) and b l B (£ e ) correspond to the ac- 
celeration of the station B and its derivative at the emis- 



sion time. The four r.h.s. terms in (24) are respectively 
zeroth, first, second and third order terms in c _1 (since 
tv B = ct (v B /c), t 2 a B — (c£) 2 (flg/c 2 ) and so on). 
Solving pi} leads to 
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(25) 
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where e = ± and with (quantities at time t e ) 
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where the third order terms in the coefficient of dt have 
been discarded since there is no zerot h o rder term in the 
coefficients of both dt e and dn k . In (28), dn k has to be 



known up to the c" 1 order only, since there is no zeroth 
nor first order term in c _1 in its coefficient. This is given 
by equation ( 26 1 , written at the relevant order 



\cn 



j] dt + [v\ - v B ] dt e + ctdn 1 = 0. 



(29) 



e = + corresponds to a photon travelling from A to B 
{tjy < t B ) while e = — corresponds to a photon travelling 
from B to A (t A > t B ). 

Considering now two emissions, at times t e and t e + dt e , 
the transfer times are t and t + dt, while the "directions" 
are n*and ri l + dn % . Differentiating t 3 b B (t e ) with respect 
to t e leads to a c -4 term, since B has at best a satellite 



like velocity. Differentiating ( 22 ) , (23), ( 24 1 and retaining 
only terms up to c~ 3 leads to 



Inserting ( 29 1 into ( 28 1 to eliminate dn leads to an equa- 



tion which can easily be solved in dt/dt e . It turns out 

(3/2) 1 

that x ph leads to terms of order c or more. Hence, 
this term is not needed. Finally 



with 





(1/2) 


(1) 


(3/2) 


dt 


dt 


dt 


dt 


dt e 




\ 

dt e 


dt e 



(30) 



cn — v 



B ~ t a B 



d^x 



ph 



dt 



-t 2 b B 

2 B 



a (3/2; 
a x 



ph 



v\ - v l B - ta l B - -t 2 b B 



v 



dt 



JX) 

k a X ph 

4 dx k A 



dt 



dtp 



ctS 



,,<• , ° x ph 

dn k 



(3/2) 1 



d x 



ph 



dn k 



dn k 



0. (26) 



All the coefficients are taken at the emission time t e - 
Since both n l and n l + dn 1 are normalized, one has also 



n k dn k = 0. 



(27) 



Solving the four equations (26), \21\ gives the four quan- 
tities dt/dte and dn k /dt e - In the time transfer problem 
we are dealing with, only dt/dt e is required. 



Making the scalar product of (26) by n l , it turns out 



(1/2) 

dt 

dt e 

(classical Doppler effect) 



A U B 



(i) 
dt 

dt e 



A B 



J Ak V B 



(3/2) 

dt 
dt„ 



ZlAlkVB 

c \ c 



: ta\ 



(31) 



(32) 



1 A 2 b\ 



v\ n k d^x 



n ■ 

k 
ph 



J Ak 



.ta k B 



ph 



dx k A 



jV l B -v A l 9x 



ct 



ct 



ph 



dn k 



(33) 



This expression has now to be inserted in (21 ). Howe ver , 
v 2 B and Wb are considered at the reception time in 



(where dt A = dt e ), while all the coefficients in (31 1 
are considered at the emission time. It turns out 



V% (t e + t) 



2tv B a B 



0(c 



(34) 



W B (t« + 1) = W B + tv B diW B + 0(c- 2 ) (35) 

all the terms being considered at time t e when not pre- 
cised. One obtains finally (suppressing the zero indexes) 



dr B 
dr A 



i U B 



1 + n 



v% , W a -Wb 



+n 



-tv 



2c 2 

i v B-K 



C 

,2 „,2 



n — -n k — 

c c 



(36) 
, cfa' n 



v z B , 



2c 2 



(3/2 

B „2 



where the dt/dt e terms of orders (1/2) and (1) are written 
explicitely to exhibit first and second order terms, cor- 
responding to classical Doppler effect, special relativistic 
contributions and Einstein's gravitational effects. All the 
c -3 terms have been grouped on the third and fourth line 



of equation ( 36 1 . Remark that OiWb in ( 35 1 and ( 36 ) cor- 
responds to a B in the case of free fall motion (i.e. if B is 
a satellite). 



2. Useful formulae 

Let us consider measurements made from a ground- 
based station (G), emitting a photon towards a satellite 
(S). This photon is reflected at (S) and come back to (G). 
To each event (emission, reflexion, return) is associated 
a proper time measured on the involved clock. Let t g be 
the proper time on (G) at emission, t$ the proper time 
on (S) at reflexion, and Tq the proper time on (G) when 
the photon comes back. Two problems may be tackled : 

- the " distance problem" , involving Tq — Tq ; 

- the "offset problem", involving [Tq + Tq) /2 — Tg. 
Obviously, the (gravitational) theory doesn't give the 



values of this last quantity, but allows to calculate dr G 
and drs as soon as dr G and the orbits of the ground- 
based and the onboard clocks are known. Reciprocally, 
the measurements of these quantities lead to constraints 
on both the orbits and the real (unperfect) clocks used 
in the time-measurement process. 

Hence, let us define the "distance problem" and "offset 
problem" quantities (observables) by, respectively 



DP = 



OP = 



d 

dr G 




d 




dT G 





TS 



Remark that since both (G) and (S) data are available, 
these quantities can be evaluated as functions of both 
Tq, Tg or Tg. However, in relevant applications, the (G) 
clock has a precision sensitively better than the (S) one. 
Hence, it is well-suited using the (G) clock as reference 
time. 



3. Explicit expressions of DP and OP 

Let us calculate DP and OP as functions of the in- 
volved quantities at the emission (e = — ) or reception 
(e = +) time t G : 



DP Oh) 

op Oh) 



dT G (t e 

dT G 



1 dr t 
2dr~, 



G 



(37) 
(38) 



In the relevant applications, the ground based station is 
the Earth. It is legitimate to consider that, during the 
time transfer, (1) the Earth is rigid and uniformly rotat- 
ing, (2) the Earth gravitational field only is acting (see 
discussion in section III I. Hence, one has, from (20 1 



dT r , 



dt; 



o 



This shows the shift between the emitted and coming 
back frequencies of the photon doesn't require an Earth 
model (for both potential and rotation). The interpre- 
tation is obvious, but this doesn't explicitely appear in 
formulae available in the Blanchet et al paper [?]. 
Let us put 



X' 
V i 
A i 
B' 
D 
iV i 

pik 



U G 



G 



- x s 

So 



(39) 



bh ~ b l s 



X l /D 

§ ik _ N i N k 



One finds, for the involved quantities 



For exhaustivity, let us point out the ratio of the ground 
emitted v G and the ground received v G frequencies is 
simply related to DP by 



dT+ 

= 1 

dTn 



DP. 



G 

dr s 



dT r 



= 1 



(*g) = 1 



,(1/2)' 

cT 1 D 



,(1/2)' 

cT 1 E 



c- 2 D 



c- 2 E 



,(3/2) 
3" 3 D 



,(3/2) 

-r 3 E 
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with 

(l/2) e 

D 

(if 
D 

(3/2) e 

D 



and 

(l/2) e 

E 

(if 
E 



(3/2) e 

E 



2N i V t 

2 (1 + e) (N l V l f - 2eV i V l - 2eDN l A 1 



N'V 1 



v k G v k G + (N k v k G ) 



+2N%P u v G V l - 2 (N l V l f N k v k s 
+D [2v%a k G + 2N k v k G N l a l G + P lk V l a G + 3V l A 1 ] 
- [(3 + 4e) N*V l + 2N l v G ] [DN l A l + P lk V l V k ] 
+D 2 N l [b G + B l ] 



2 
D 



k 2^ 

P lk V l c 2 ^ ph -2N l V l N kC ph 



cdt 



.1 N k pit yi^ph _ 2N k v l ^ dX ph 

D dn l s dx l a 



N l V l 
T 



+ e) V l V l + (l + e) (N l V l ) 



DN l a G 



+eDN l a s + W G -W S 



1 



N i v%V k v k q + -N l V 



(N k v k ) 



-2(1 +e) N l V l P kl V k V l 
1 



+D 



+D 



yV s + (l + e)^ (2a 4 G -4) 



N i V t N k a k s + N l v s N K a s 



i n A \rk „k 



k \k 



-2(1 +e) DN l V l N k A 



+D Z 



-l N % + (l + e)N% 



+N i V i (W G - W s ) + eDv s d,W s 



1 

D 



P lk V l c 2( x ph - N l V l N k 



2«« 
CO X 



cdt 



D dn l b dx l s 

The quantities D, N,v Gl vs,a Gl as,b G ,bs are taken at 

the time t G , while the photon first order term ^x p h 
and its derivatives have to be computed for t — +D/c, 
n = +Ax/D and x$ (at time t G ) for both e = ±. Let us 

9 (i) 

recall that c x p h and ct are zeroth order quantities. 
From these expressions, it turns out that 

(1) e (3/2) c 

OP(t e G ) = c~ 2 OP +c- 3 OP 



with 

(i) 6 
OP 



(3/2)' 

OP 



+ DiVV G - V^ G + W s 



—N i V i (N k v h s ) 2 - ^V' l V l N k V k + (JVY') 3 



+D 



v l G a G - eV (A 1 + a G ) + -N l V l N k a k s 



CONCLUSION AND DISCUSSION 



In this paper, one has proposed formulae involved 
in two way time transfer problems in a form directly 
adapted for practical uses. All these quantities are ex- 
pressed as functions of quantities defined at a unique 
(emission or reception) ground-based time. Considering 
foreseeable increase of clock's accuracy, the formulae are 
given including the J2 contribution, but can be restricted 
to the spherical case making J 2 = 0, which is sufficient 
for the time being applications. 

In this paper, the chosen unique time is a ground time 
because ground based clocks have generally precisions 
sensitively better than on board clocks. However, in the 
case where one needs to make some time transfer between 
two different ground based stations, distant in such a way 
there is no spacecraft in common view, the link will in- 
volve an on board clock, with the resulting limitation in 
precision. Let t\ and ti be the times of the links between 
the spacecraft and the first and second ground stations. 
Since the spacecraft motion is considered during the time 
interval [ii,^], the crossed cdtdx 1 term in S is a c" 4 
order term. Hence, the onboard proper time formula re- 
duces to 



Ut 



2w 

^2 



+ 0(c- 4 ) 



leading to the spacecraft proper time interval 



At (t!,t 2 ) = t 2 - h + 1 ^ ( W ~y) 



dt + 0(c- 4 ). 



Let us enforce the fact that, since the time interval [t±, £2] 
is finite, the orbit used in the previous r.h.s. member inte- 
gral cannot be developped around the emission time, and 
should even include all the perturbations acting on the 
satellite's motion. On the other hand, a classical descrip- 
tion of this orbit is sufficient, since relativistic corrections 
generate c -4 contributions to At. 
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